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Abstract 

We examine the stability of charged Lovelock black hole solutions under vector type and scalar 
O I type perturbations. We find the suitable master variables for the stability analysis; the equations 

C/3 ' for these variables are the Schrbdinger type equations with two components and these Schrodinger 

operators are symmetric. By these master equations, we show that charged Lovelock Black holes 
are stable under vector type perturbations. For scalar type perturbations, we show the criteria 
for the instability and check these numerically. In our previous paper, we have shown that nearly 



5-H 

extremal black holes have the instability under tensor type perturbations. In this paper, we find 
that black holes with small charge have the instability under scalar type perturbations even if they 
have relatively large mass. 
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I. INTRODUCTION 



The braneworld scenario with large extra dimensions predicts that higher dimensional 
black holes might be produced at colliders Therefore, higher dimensional black holes 
become attracting subjects and some aspects of these have been inspected so far. For 
example, exact solutions are investigated in higher dimensions. In higher dimensions, besides 
Schwarzschild black hole, Reissner- Nordstrom black holes Q and rotating black holes jsj, 
various solutions are found: black ring solution 14 1, black di-rine 

so on. For these solutions, from the standpoint of black hole creations, it is important to 
examine the stability of such solutions because stationary solutions with the instability are 
not attractors of time evolution. This suggests that such black holes should not be realized. 

So far, various stability analyses for black hole solutions have been performed. One of 
the most notable analysis is that of Tangherlini-Schwarzschild solutions by Kodama and 
Ishibashi [7|. They have derived master equations for all type perturbations. These are 
Shrodinger type equations and they have shown that these Shrodinger operators are all 
positive definite using the S-deformation approach which they have developed by Friedrichs 
extension. These results show that Schwarzschild black holes are also stable in higher di- 
mensions. They have also examine the stability of higher dimensional Reissner-Nordstrom 
black holes [8|. By the S-deformation, they have also shown that this charged solution is 
stable under tensor and vector type perturbations. For scalar type perturbations, it has 
been shown that this black hole is stable in 4 and 5 dimensions. For this solution, the stabil- 
ity has also studied numerically and it has been found that black holes with large negative 
cosmological constants and large charge are unstable in more than 7-dimensions [9|. On 
Myers-Perry black hole solutions, in D > 6, there found the instability for singly rotating 
solutions when the spin parameter is large enough [ji]] . In even dimensions, stability of near 
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horizon geometry of rotating black hole with equal angular momenta are investigated 
It has been suggested that scalar mode for base space has the instability. Recently, the 
stability of black ring solution is examined by using local Penrose inequality and it is shown 
that the fat branch is unstable [l^ . 

The stability analyses we have introduced above are all premised on Einstein theory. In 
fact, the stability of black hole solutions have been examined mainly in Einstein theory. 
It is as important as such analyses to investigate the stability in more general theories. 
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In 4- dimensions, Einstein theory is characterized by two properties; the action has the 
general coordinate covariance and equation of motion consists of metric, the first derivative 
of metric and the second derivative of metric [kJ. Then it is natural to extend the four 
dimensional gravitational theory to higher dimensional one keeping these two properties. In 
higher dimensions, the most general theory which satisfies above two features is not Einstein 
theory; it is Lovelock theory [l^j. Then it is important to generalize the stability analysis 
of black hole solutions in Einstein theory to these in Lovelock theory. 

Same as Einstein theory, a spherical symmetric solution is known in Lovelock theory 
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16j | . This solution is called as Lovelock black hole solution. For this Lovelock black hole 



solution, the stability has been analyzed in 
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19[. In these papers, it has been shown that 



black holes with sufficiently small mass are unstable under scalar type perturbations in odd 
dimensions and unstable under tensor type perturbations in even dimensions. This critical 
mass differs with dimensions and Lovelock couplings. These instabilities become stronger 
as wavelength becomes smaller, and the time scale of the instability converge to in small 
scale limit. Under vector perturbations, this solution is stable in all dimensions, which is 
independent of mass. 

Since black hole creations originate from protons at colliders, it is also important to take 



account o 
is known 



Maxwell-charge. In Lovelock theory with U (1) field, a charged black hole solution 



15J ; this has spherical symmetry and a time-like Killing vector, and this Reissner- 



Nordstrom like solution is called as charged Lovelock black hole solution. Then, in this 
paper, we'd like to extend the stability analysis for Lovelock black hole solutions to charged 
Lovelock black hole solutions. For this charged solution, stability analysis under tensor type 
perturbations has been examined by us and we have shown that black holes are unstable if 



they have nearly extremal mass 



20| | . In this paper, we extend our previous discussion to 



vector type perturbations and scalar type perturbations; that is, we derive master equations 
for these type perturbations and examine the stability using master equations. 

The organization of this paper is as follows. In section HI1 we review Lovelock theory, 
present the charged Lovelock black hole solutions and check the behavior of these solutions. 
We mainly concentrate on asymptotic flat branch. In section IHH we review the analysis for 
tensor type perturbations 20|. We examine tensor perturbations and show the criteria for 
stability under this type perturbations. In section IIVt we derive master equation for vector 
type perturbations and show that there is no instability under this type perturbation. In 
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section |V] we concentrate on scalar type perturbations. We show that master equations can 
be summarized as a Schrodinger type equation with two components, and using this equation 
we present criteria for stability. In section IVI[ we numerically examine the conditions for the 
instability presented in section II III and [V] In this paper, we only check in 5 — 8 dimensions. 
In the final section IVII[ we summarize this paper. 



II. CHARGED LOVELOCK BLACK HOLES 



In this section, we introduce Lovelock theory and present charged black hole solutions 
in Lovelock-Maxwell theory. These solutions are expressed as the roots of the polynomial 
equation and we confirm that one of the roots is asymptotic flat. For this asymptotically 
flat root, we briefly check the behavior, singularity and horizons. 



A. Lovelock-Maxwell System 



f. Q, 



In Ref. [141 ]. D. Lovelock have constructed the gravitational theory whose equation of 
motion consists of the metric, the first derivative of the metric and the second derivative of 
the metric. The Lagrangian for this theory is 

£ Lovelock = _ 2 A + fiiR 

h 2 m ™ iC; 2 (n - p) Ki pi <v K ^ K ^« m 

where A corresponds to a cosmological constant and /3 m s are arbitrary constants which we 
call Lovelock couplings. We add the coefficients (2m)!/2 m m np=i~ 2 ( n — P) f° r convenience. 
In the above Lagrangian, n is related to dimension D as n = D — 2 and k corresponds to 
the maximum order defined as k = [(D — l)/2] where [x] is the Gauss symbol. There exists 
maximum order k due to the antisymmetric property of Sl^S^d^ ■ ■ ■ 5p™ . When we fix the 
maximum order k, the dimension is restricted as n — D — 2 — 2k — 1, 2k; for example, the 
second order Lovelock theory is the most general in n = 3 or n = 4, and the third order one 
is in n = 5 or n = 6. By the ambiguity of the overall factor of the action, we take the unit 
0x = 1 in this paper. 

In this paper, we want to concentrate on Lovelock-Maxwell system. This system is 
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described by the action 

S = J d D x^gC Lovelock - j d D xV=g~F^F>*», (l) 

where is the field strength of Maxwell field A^. In the action (0Q), the dynamical variables 
are 9 [iv an d A^. The variations by these variables lead 

Q» = T» , (2) 
F" v , v = , (3) 

where Q^, which we call Lovelock tensor, and T^, which means energy momentum tensor 
for U(l) field, are defined as 

g» = as; + v - Irss 

_ (2m + 1)! 0rn A^A A i A^i A A 2 . . . X<x m ] d ...p «mPm (A) 

^2 2 m U P =i {n-p) 
T» = F fl \F uX - \F Xp F^5; . (5) 

The field strength is defined as F = dA, then F pv must satisfy the identity 

dF = F [fiu . x] = 0. (6) 

The above equations (J2J), (EJ) and (jSJ) are our basic equations. 



B. Charged Lovelock Black Holes 



15|. We 



For the basic equations, black hole solutions with two parameters are known 
assume the static spherical symmetric metric with spherical symmetric electric field 

ds 2 = -f(r)dt 2 + l/f{r)dr 2 + r 2 llj dx i dx j , (7) 
F tr = E(r), other components = . (8) 

In these, 7^ corresponds to the metric for S n . 

We can easily check that these ansatz satisfies fl6]). Then we concentrate on the others 
and these lead following equations; 

= d r {r n E{r)) = , 

Gt* = Tt, QS = T/ =► (r" +1 P[^])' = , (9) 
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and the other components are identical. In ([9]), ip is related to f(r) as / = 1 — r 2 ip and V[ip] 
is defined as 



m = E 



m=2 



/3 



m 



2A 



n(n + 1) 



(10) 



The second equation of is derived by a derivative of the third equation with the first 
equation, so we only consider the first and third equations. The first equation can be easily 
integrated and the result is 

E(r) = y/n(n - l)Q/r n . (11) 

In this equation, ^/n(n — 1)Q is an integral constant and this constant corresponds to the 
charge, which can be seen from the behavior of E(r). Substituting (ITTj) into the third 



equation of 



), we can gain (r n+1 V 



V\ 



M 

r n+l 



[n — l)Q 2 /r n , or integrating both sides reads 
Q 2 



r 



2n 



M(r) 



(12) 



where M. is an integral constant. We will see that M. corresponds to mass when checking 
the asymptotic behavior of the solution 2l|, |22|. 

We must solve the polynomial equation ( fl2l) for solution of Lovelock-Maxwell system. In 
order to solve the polynomial equation (1121) . in this paper, we assume some conditions for 
Lovelock couplings /3 m and Ai for simplicity. First we consider mass of black hole is positive, 
that is, Ai > 0. Second we set cosmological constant A = 0. In A = 0, as we will see later, 
there must exist an asymptotic flat branch. Third, for simplicity, we assume the positivity 
of Lovelock couplings, that is, 

Pm > (m > 2) . (13) 



C. Asymptotic Flat Branch 

Because ([12]) is k-th order polynomial, the polynomial equation ([12]) should have at most 
k solutions. However, assuming all Lovelock couplings are positive and A = 0, one of the 
roots corresponds to an asymptotic flat solution. 

For instance, we see the above statement when k — 2. In this case, (Tl2|) reduces into the 
2nd order polynomial equation, so this can be easily solved as 

-1 + + 2foM/r^ - 2/3 2 Q 2 /r 2 ™)//3 2 

mr)={ , . (14) 

(-1 - v 7 ! + 2^M/^+ 1 - 2f3 2 Q 2 /r 2n )/(3 2 
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FIG. 1: We introduce the graphical method for finding the roots of (|12p in this figure. The solid 
curve corresponds to y = V[ip] and the dotted horizontal line is y = M(r) with fixed r. The cross 
points in this figure are the solutions of the polynomial equation (|12p for this r. If we want to 
consider the roots for other radii, we draw the corresponding dotted line y = M(r) and see the 
cross points. 

Let's consider the limit r — > oo. The first root behaves as ip — > M./r n+1 and the second 
converges as ip — > — 2//3 2 . Therefore, the function f(r) = 1 — r 2 tp(r) behaves as / = 
1 — Ai/r n ~ 1 for the first root and f = 1 + 2r 2 //?2 for the second one. This shows that 
the first branch expresses an asymptotic flat solution and the other is an asymptotic AdS 
solution. 

Regrettably, we can not write the roots of ([121) with general k explicitly. Nevertheless, 
we can understand the existence of an asymptotic flat solution as the solution of (I12p . In 
order to check this statement, we want to introduce the graphical method. In Fig{TJ we 
line y = M(r) with fixed r and y = V[ip] in ip — y diagram. The former is a horizontal 
line because M depends only on r and we fix r. The cross points in FigfTJ are roots of the 
polynomial equation ( TT2|) for this fixed r, and if we want to find the roots for other radii, 
we move the horizontal line following the value of M(r) and check the cross points. 

For this method, it is important to check the behavior of M(r) = Ai/r n+1 — Q 2 /r 2n . 




Then, like Figj2j M(r) becomes 



at r = r = (Q 2 /A^) n_1 , takes a maximum value M, 



max 
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FIG. 2: We plot M(r) in this figure with n = 3, M = 1 and Q = 1. ro and r max is defined as 
M(ro) = and M'(r max ) = respectively. M max corresponds M(r max ). We can recognize the 
manners of roots of the polynomial equation (|12|) by this behavior and Fig[TJ 

r ma x= (2n/(n + l))"^ 1 r an d behaves as M ~ Ai/r n+1 in the asymptotic region. 

As mentioned above, M(r) is positive in r > r . While M(r) is positive, from FigJTJ the 
polynomial equation (fT2|) has only one positive root because V[ip] satisfies V[0] = and is 
a monotonically increasing function in ip > under our assumptions ( TT3|) . This positive 
root expresses an asymptotic flat solution. To confirm this, let us consider the behavior of 
this ip when r — > oo. Because M(r) converges to like Ai/r n+l as r — > oo, the positive 
cross point in FigJTJalso converges to 0. In detail, eq.f lT2|) with if; ~ shows that this root 
converges like ip ~ Ai/r n+1 . Then, for this branch, f(r) behaves as /(r) ~ 1 — M./r n ~ l in 
the asymptotic region. The metric ansatz ([7]) with this asymptotic behavior shows that this 
positive ip corresponds to an asymptotic flat solution. This asymptotic behavior of f(r) also 
explains that Ai corresponds to ADM mass. 

In the last of this subsection, we briefly check the behavior of our asymptotic flat root 
not in the asymptotic region. Let us consider by the graphical method with Fig{T] again. 
Because M(r) behaves as — > M max — > when r moves oo — > r max — > r , our ip{r) varies as 
— > ipmax 0. When r becomes smaller than r , M(r) becomes negative and so ip(r) also 
takes negative values. 
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D. Singularities 



In if) < or in r < ro, our asymptotic flat solution has the curvature singularity. To 
confirm this, we examine the Kretschmann invariant 



This value diverges at r = 0; and the singularity also exists where /' or /" diverge. For 
example, /' has the term like r 2 ip' . Form the the derivative of (fl2l) . this can be estimated 
as r 2 if)' = r 2 M' '/d^V '. Then, besides r = 0, R^ v \ p R^ uXp also diverges where the derivative 
of V[if>] with respect to if) becomes 0. If V[if)} takes a extreme value at ip , because V[if)] is 
monotonically increasing in if) > 0, such ip must be negative; that is, there is the singularity 
at r s (< ro). If V[if)] is monotonically incresing for all if), there is the curvature singularity at 
r = 0. Therefore, whether V[if>] has extreme values or not, our asymptotic flat branch has 
the curvature singularity somewhere in < r < ro. 

E. Horizons 

Singularities must be wrapped by the event horizon from the standpoint of cosmic cen- 
sorship. In this subsection, we consider horizons and present the condition for existence of 
horizons. 

Our asymptotic flat solution has the event horizon at f(r) = 0. This branch also satisfies 
V[if)\ = M(r), so horizons can be determined from 



where r# are horizon radii and ipn is defined as ipu = ip(rjj). The first equation shows, if 
horizons exist, the corresponding ipn must be positive. Our if)(r) is positive in r > ro, then 
thi if exists, must satisfy rn > tq. As we have emphasized, the singularity exists somewhere 
in r < r . So we do not worry about the naked singularity if (fT5j) has roots. 

Here, we consider the criteria for the existence of roots of eq. (1151) . Eliminating ipH by the 



R.uXpR^ = f" + In— + 2n(n - 1) 



(I-/) 2 



= 1- rjjif) H 
V[if) H ] = M{r H ) 



_ M 




(15) 
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FIG. 3: The solid curve corresponds to y = a(r) and dotted line is y = M.. In this figure, y = a(r) 
takes an extremal minimum at r = r ex . Because of f 1 1 6 [) . the cross points mean horizon radii m- 

first equation of f fT5|) . the second equation becomes equation for r# as follows; 

M Q 2 



V [l/rl] 



„n+l r 2n 
H ' H 



^M = ^ + (V + £ = a(r H ) . (16) 

H \ m= 2 J 

The first term of a(r) is negative power of r and its coefficient is positive. Under our 
assumption ffTS]) . because n = 2k or 2k — 1, the other terms are positive power of r and 
their coefficients are positive. Therefore, y = a(r) behaves as FigJ2l ot{r) diverge near r = 0, 
takes extremal minimum at r = r ex and monotonically increase in r > r ex . Then, we can 
denote that f|T6|) has two roots when A4. is larger than Ai ex where 

M rI = a(r ra ) = + 4- Y, ^:»-A . (17) 

ex \ m=2 / 

In the two roots, the larger one corresponds to the outer horizon and we call this r out 
hereafter. Note that a(r) only depends on Q except for Lovelock couplings. Therefore, r ex 
is determined when we fix Q; this shows that M. ex depends only on charge. 

Finally, we check the behaviors of some functions for the later discussions. First, we 
examine the behavior of M(r) outside of r out . The l.h.s of the first equation of (TL6l) is 
monotonically decreasing function and the r.h.s is monotonically increasing in r < r max and 
monotonically decreasing in r > r max . Then, when f|T6l) has two roots, it is forbidden that 
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both of them are smaller than r max ; at least, the larger root r out must satisfy r out > r max . 
Hence, 

M'(r) < (r > r out ) . (18) 

Then, in FigJTJ while r > r out , the dotted line y = M(r) falls monotonically as r becomes 
larger and so ip(r) decreases monotonically in this region. Therefore, from the relation 
i/j'd^V = M' reads 

djVWl > (19) 

when we consider the outside of r mt . 



III. TENSOR TYPE PERTURBATIONS 

Thanks to the spherical symmetry of background (JTj), tensor-type, vector- type and scalar- 
type perturbations are decomposed and we can examine them separately. We have already 
examined the tensor perturbations in [^(J. In this section, we review our previous analysis 
briefly. Note that we only consider the case when there exist horizons; that is, M. > M. ex . 



A. Master Equation 



Under tensor type perturbation, there is no perturbation for Maxwell field. Then we only 
consider the gravitational perturbations 





' \ 


\ 












(20) 




v o 


r 2 07-, J 





In this expressions, <fi corresponds to the master variable. Hj is the tensor harmonics which 
is characterized by the traceless condition T\ = 0, the divergence-free condition = 
and the eigenequation 7~i/ k \k = —(£(£ + n — 1) —2)%j. Note that | is the covariant derivative 
for and i is the integer which satisfies i > 2 . 



Using above metric perturbations, the first order equation 5Qi 



T'<f) - fT'cf)" - f 2 T 



Ar 2 fT>) 
r 2 fV 



+ n-l)f_ T „ 



(n — 2)r 



leads 
= , 



19] 



(21) 
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where T(r) is 

T(r)=r n - l d 1 pV[ij] J (22) 



which is always positive in r > r out due to (U9I) . 



For this equation, as we have shown in 



19|, 



20|, there exist ghost like instabilities if T' has 



negative regions. For example, the coefficient of kinetic term in ()2ip is proportional to T', 
so this term has the wrong sing while T' is negative. Then, here we also assume T'{r) > 
in r > r out for avoiding the ghost instability. 

Under the ghost-free condition, we can change the normalization of <fi as \&(r) = 
4>{r)r\jT'{r). Using this variable, converting r to r* which defined as dr* / dr = 1// and 
Fourier transforming like \& — > ^e luJt , (12ip is recast as 

m> = ou 2 ty , (23) 

where 

% = -d 2 r , + V g {r) 

+ n-l)fT" , 1 a2 



W= \n-2)r T + 7VF d ^ Vr) - ^ 
Eq. (l23p is a Schrodinger type equation and its eigenvalue is u 2 . Then, if this Schrodinger 
operator % has negative spectra, we can say charged Lovelock black holes are unstable under 
tensor type perturbations. 



B. Stability Analysis 

In this subsection, we show "there exist negative spectra if T" has negative region in 
r > r out " when T' is always positive. To show this, we define the inner product as 

/oo 
X*iX2dr* , (25) 
-oo 

and use the inequality 

(X,^X)>^ 2 -(X,X), (26) 

where u 2 is the lower bound of spectra and \ is an arbitrary smooth function with compact 
supports. From this inequality, we can show that there exist negative spectra if we find a 
function \ suc h that (\, Hx) becomes negative under our assumptions. 

12 



We assume T" has negative regions and define / as a closed set on such regions. Under 
these, we chose Xo as a smooth function which has a compact support on /. For this Xo, 
(Xoj^Xo) is evaluated as 

(Xo,^Xo) = fdr*[-xld 2 r , X o + V g \xo\ 2 } 

= I \d r *Xo ~ f± In (rVr) X o\ 2 dr* + £(£ + n - 1) jf j^y^dr* . (27) 

In this calculation, we use Gauss divergence theorem and neglect boundary terms because 
Xo is smoothly connecting to at dl. In ( 1271) . the first term must be positive and the second 
integral is negative because we assume T' > in r > r out and T" < on /. Therefore, taking 
£ — > oo, (xo? Hxo) must become negative. Because of (}2T)j) . this means negative spectra exist 
in sufficiently large £ modes. Then we can declare that black holes are unstable if T" takes 
negative values somewhere in r > r out . 

Inversely, it can be also shown that charged Lovelock black holes are stable if T" is always 
positive in r > r out . As shown in Ref. it is sufficient for the stability to show that ($, 
is positive for v $ e C^°(r*). We can check this criterion by the same calculation of ( 1271) and 
the positivity of T", so we can say black holes are stable if T" is always positive. 

We want to summarize this section. For avoiding the ghost instability, we must assume 
T' is always positive. Under this assumption, charged Lovelock black holes are stable if and 
only if T" always takes positive values in r > r out . Hence, what we have to do is a probe of 
the behaviors of T' and T", and we will check these in the section IVI1 

In the end of this section, we want to comment about T' in the 2nd order Lovelock 
theory. We have already shown in our previous paper [2^], there is no ghost instability for 
this case. This can be checked by the direct calculations; eq. (|T4l) leads T(r) = r n ~ l d^V[^)\ = 
r n_1 ^/l + 2/3 2 M(r), then T is calculated as 



T 



'I + 2/3 2 M(r 
IV. VECTOR TYPE PERTURBATIONS 



( r2 _l) + ( r2 _3)__ + 2-^- 



> . (28) 



In this section, we examine the vector type perturbations. We here also assume A4 > A4. ex 
for horizons and X" > for no ghosts under tensor type perturbations. We only consider 
the perturbations in r > r out . Under these assumptions, we derive master equations for 
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vector type perturbations. Using this equations, we show that charged Lovelock black hole 
solutions are stable for vector type perturbations when T' is always positive. 



A. Gravitational Perturbations 



Firstly, we'd like to consider the metric perturbations. In this paper, we use the Regge- 
Wheeler gauge in which metric perturbations are expressed as 



/ 













\ 







(29) 
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y sym sym 

In these, Vj is the vector harmonics which is characterized by the transverse condition = 
and the eigenequation V^ l \i = —K, v Vi with k v = £(£ + n — 1) — 1 (£ > 1). 

For vector type perturbations, other than 5Qt-, $G r l and 5Q^ are trivial. From the above 
metric, we can calculate the non-trivial components as 



SOS 

sg r i 



n-l))T> f 
rii — 



2{n - l)r n + 2 

K v -(n- 1))T 
2(n - l)r n + 2 



ho- 



2 r n+2 

1 r 3 T 



E 1 



V 



1 



-jhi + (fT%)' 



2(n - 1 

This is almost same as the results of 



2r n+2 f 



+ —h 2 



(30) 



19[ except for the background electric field. This gap 



mainly arises from the difference of the identity 

{{rf + 2{l-f))T}' = 2 K —L^ 0) 



B. Perturbation of Maxwell field 



Next we examine the vector perturbations of Maxwell field. We start from the perturba- 
tion of vector potential 

5A„ = (0, 0, CVif , (31) 

where V, is the vector harmonics. Note that C is gauge invariant under U(l) gauge because 
there is no gauge freedom for vector perturbations. Using above SA^, we can easily calculate 
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the first order of the field strength as 



5F ti = CVi, 5F ri = C'Vi, 5F tj = C (% - %) , otherwise = 



(32) 



Here, we derive the evolution equation for first order variable C from Maxwell equations. 
It is easy to check the above field strength satisfies the identity SF^ U .\] = 0. Therefore, 
5(F^ IU ' U ) = is important for the evolution equation. In these equations, \i = t,r components 
are trivial and /i = i components read 



jC - ±fi (r-MC) - r*E ( (£)' - M + ^±^>C = . 



(33) 



In order to gain evolution equations for gravitational field, we must calculate the first 
order of the energy momentum tensor (jSJ). This tensor consists of 5F^ U , Sg^ and background 



variables, so this tensor can be calculated from 
as follows; 



(J7|) and (jSJ), and the results are 



5T r l 



r E 2 



V 
V 



other components = 



(34) 



Then, from ( 130]) and ( )34|) . the first order Lovelock equation bQ^ = 5T^ U reads 



V 



2(n - l)r™+ 2 
V 

2(n - l)r n+2 



(k v - (n - 1)) h\ 



(k v -{n- 1)) h 2 



f 



ht\' h 



>- 



1 x\M(h^-h 



2r n+2 j 



-jk + (fT'h 2 )' = 



fr 2 



(35) 



C. Master Equations 



Now we are position to derive master equations from (1351) and (|33|) . First of all, we treat 
the third equation of fl35|) . From this equation we can define a new variable <fi as 

1 • 



(36) 
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Substituting (136]) into the first equation of (135]) and integrating this with respect to r reads 
W + Cl(1) .^_i),_ r ^(^_^ . ,37) 

Here we use (r n E) = const, and C\(t) is a constant of integral. Same as this, substituting 
( 136]) into the second equation of ( 135]) and integrating this with respect to t reads 



2 (En C + 0,(r) = - ,T ( Y - £ ] . (38) 



Comparison ( 137]) with ([38]) shows that Ci(t) = C*2(r) = const, and this constant can be 
absorbed into <fi. Therefore the three equations (|35[) are reduced into one equation 

W = =^=f^ - r*T ( (-^) ' - ^) . (39) 
Same as this substitution, ( 136]) makes Maxwell equation ( 133]) 

yC - - ^5* (r- 2 fd,.c) 

+ 4 ( 2(E ,.- )g - — - - % ) + ^L±ipl) g = o . (40) 
ri \ n — 1 / 

As we have seen, determine the perturbation of gravitational field and C does Maxwell 
field. Therefore, these are the master variables and (139]) and ( jjOj) are the master equations 
for vector type perturbations. 

Finally, we'd like to alter these two equations into a Scrodinger equation with two com- 
ponents. To do so, we must change three points: Firstly, we change the normalization 
as 



rVT y^-(n-l) 1 ; 

In this, we use the assumption that T' is always positive. Secondly, we switch radial coordi- 
nate r to r*. Finally, we Fourier transforme like \l/ — > tye luJt and £ — > £e la; *. Then, ([39]) and 
10]) become a Scrodinger equation with two components as 

(42) 

\ C J \ C J 

where 





« = W ^ | ,43) 
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and 



\r < \ Kv-(n-l) fT' r— 2 1 



VJr) 




r\JT> 

n-l) fVr 



2 {n - 1) r ("+ 2 )/ 2 T ' 



V em {r) = fa + (n- + 2{Er n f-l If + r -^ 2 d 2 ,r^ 2 . (44) 

Note that these equations are not decomposed due to higher curvature collections. In Ein- 
stein limit, owing to T(r) = r n_1 , the above potential matrix can be diagonalized by con- 
stant eigenvectors. This indicates that our Schrodinger equation can be decomposed into 
two equations by taking suitable linear combinations of ip and (. In fact, we can do so by 
using the combination a±ip — ( with 



- 1)M ± vV 2 - l) 2 M 2 + 8 fa - (n - l))(r n E)' 



a± = - — -. — — , (45) 

2(r n E)y/2( Kv -(n-l)) 

which is consistent with Against this, because T(r) is more complicated in general 
Lovelock theory, we must consider the above coupling system. 



D. Stability Analysis 



In this subsection, we show that the Schrodinger equation (1421) has no negative eigenvalue 
states. 

In order to show this, for \I/ = (\&, () T , we define the inner product as 



dr* [*;w a + crc 2 



(46) 



Here, we prove that charged Lovelock black holes are stable for vector type perturbations. 
In order to show this, we prove % is an essentially positive-definite self-adjoint operator. For 
this, because % with C^°(r*) x C^°(r*) is a symmetric operator, it is sufficient to check this 
operator is positive-definite Q. 

We assume \l/o £ Co°(' r *) * C^°(r*), then (^/q,'H^q) can be estimated as 



dr* 



*X2 



+ V g \^ \ 2 + VMCo + ^oCo) + ^emlCol 
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dr* 



a c n " 2f c 

°V*V>0 7; S>0 



2 r 



~ 1} ^o| 2 + 2(Er") 2 ^-|Co| 2 

n — 1 rl r n+ij 



+(«„ + (n- l))^|Ci 



2(«« - (n- 1)) n.h/T', ,*>. , , 



dr* 



n — 1 v ' y r ( n+2 )/ 2 T 
2 



9 r ^ + (5* ln(rVT'))Vo 
V2(Er n 



Or*^0 ^ C,0 



+ 



rT 



„n/2 



2 r 
2 



. k„ - [n - 1) 
Co - \ : vT'^o 



+ (^, + (n-l))4lCo| 2 



n — 1 

>0 , 



(47) 



where we use Gauss theorem in the second equality and neglect boundary terms because 
^0 and Co are in C£°(r*). This calculation shows that % with C£°( r*) x C^°(r*) is positive- 
definite. Then, since % with x is essentially self-adjoint, % can be uniquely extended 
to a positive-definite self-adjoint operator, so there is no instability under vector type per- 
turbations. 



V. SCALAR TYPE PERTURBATIONS 

In this section, we derive the master equations and present conditions for the instability 
under scalar type perturbations. In this section, we assume Ai > Ai ex and also assume that 
T' is always positive outside r out . 

A. Gravitational Perturbations 

Firstly, we consider the metric perturbations. In this paper, we take the Zerilli gauge in 
which metric perturbations are described as 





' fH y 







\ 




sym 


Hy/f 









y sym 


sym 


r 2 Ky^ 


J 
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where 3^ is the scalar harmonics which is characterized by the eigenequation y^ 1 = —K s y 
with k s = £(£ + n - 1) and t = 0, 1, 2, • • ■ . 

For deriving the master equation, it is sufficient that we calculate bQ\ {i ^ j), 6Q%, 6Qj, 



bQ\ and 5G r r \m 



By the above metric perturbations, we can derive the following results 



19|; 



— — l -—(T'H - T'H - rT"K)y lj 



\i , 



fT 

2r n+l 
1 

2r n+l 



-K s Hx + n < r 



v 



K + r 2 K' - rH 



y 



{—k s T - nr(fT)'} H - nrfTH' + (n - K s )rT'K 
nr 2 f'T 



5Gr l 



2r n+1 
1 

2r n+1 



+ nf(r 2 T)' } K' + nr 2 fTK h 



f'T T 



y 



) H + TH' Q - jH ± 



nr 2 T 



2/ r 



T 



nr 2 f'T\ 

'^~J^ K + {n- K s )rT'K + ( nr 2 fT 4 



-nr(fT)'H + k s TH - nrfTH' 



y . 



(49) 



These are same as our previous calculation for neutral black holes except for the detail 
expression of f(r). 



B. Scalar Perturbations for Maxwell Field 

Next, we examine the scalar perturbations for Maxwell field. We start form perturbations 
of the field strength which has the U(l) gauge invariance. We describe this as 





( o xy 






SF fa/ = 


-xy o 


zy v 


(50) 




^ anti sym. 


J 





where y is the scalar harmonics. Note that 5Fij = because we can not construct antisym- 
metric tensors from scalar functions. 

The start point (150]) enables us to calculate the identity SF^ U . X ] = 0, Maxwell equations 
5{F^ V ' - v ) = and the energy momentum tensor ST^. Firstly, we check the identity. It 
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is easy to show that the components other than (/i, v, A) = (t, r, i) are trivial and the 
non-trivial component reads 



X = Y' - Z . (51) 

Secondly, we calculate Maxwell equations 6(F IJ ' 1/ . V ) = 0. From \i — t component, we can 
gain a equation 

I n—2 

-(r n E)d r (H -H + nK)- d r {r n X) + ^— Y = . (52) 
* J 

Same as this, \x — r component reads 

^(r n E)d t (H -H + nK) - d t (r n X) + K s r n ~ 2 fZ = (53) 
and fx = i components are 

d r (r n - 2 fZ)-d t (-Y) = 0. (54) 

Before calculating ST^, let us reduce the four equations (15T|) ~ (!54|) into one equation. 
From (EH), we can define a new variable B as 



B, Y = -L.B'. (55) 



r n-2 j ' r n-2 

Substituting 055p into 052 1) and integrating this with respect to r, becomes 

^(r n E)(H -H + nK) - r n X + k s B + d(t) = , 

where C\{t) is a constant of integral. We use r n E = const, in this integral. Same as this, 
( 1531) becomes 

^{r n E){H -H + nK)- r n X + k s B + C 2 (r) = . 

Here ^(r) is also a constant of integral. Then, a comparison with above two equations 
reads C\(t) = Ci{r) = const, and so this term can be absorbed into B. Therefore, 

- (r n E) (H -H + nK) - r n X + k s B = . (56) 

Eq. (1551) and Eq.( l56l) suggest that B (and gravitational perturbations) determines the per- 
turbations of Maxwell field, so we can say B is the master variable for Maxwell field. The 
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evolution equation for the master variable B can be derived eliminating X in (156]) by (l5Tj) 
and fl55|) . The result is 



d 2 B -r n d r - 
f r 



f 



-0 r B + 



(r n £) 



(#o - H + nit) = . 



(57) 



Finally, for the first order perturbations of Lovelock equations, we calculate the first 
order of the energy momentum tensor (j3J). Eq.((5]), the background electric filed E(r), metric 
perturbations (1481) and perturbations of Maxwell field (150]) yield 



ST/ = 8T r 



EX + -E\H - H ) 



y 



K 



n 



— B + -EK 

r n 2 



y 



0, 5T r 



E 



:Yy 



E 



iy 2 ^ ly I 



(58) 



Note that we use the relation (1561) in (t,t), (r, r) components and also use the relation ([55 
in (r, i) components. 



C. Master Equations 

From the first order Lovelock tensor ( I49j) and that of energy momentum tensor (15 8 j) . the 
components we concentrate on are 

T'H - T'H = rT"K , (59) 
-k s H x + n | (r - T -^j K + r 2 K' - t h\ = , (60) 

(n - K s )rT'K + (V^fL + n f(r 2 T)'\ K' + nr 2 fTK" 



-(k 8 T + nr(fT)')H - nrfTH' = 2r n+1 E {^B + ^EK^j , (61) 

( r + 27 J H ~ ' + v V " 7 ) H ° + TH '° ~ 7 Hl = 2ErB ' ' (62) 

/ nr 2 f'T\ 
(n - K s )rT'K + I nr 2 fT' + — - — J K' - nr(fT)'H 

+k s TH q - nrfTH' + 2nrTH l - °^-k = 2r n+1 E f^B + ^nEKJ . (63) 

From now on, we'd like to construct the master equation. In order to derive, we must 
define the master variable <fi and denote the gravitational perturbations by this 0. Same as 
the analysis for neutral black holes [l^, we define the master variable <fi as 

H^jty + K) . (64) 
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Then, from (159]) . (!60|) and (ISTj) . we can express ifo ; H and as follows 19| ; 
H 



rT" 
H + —K 



nf 2nj 



inrfE 2 
.4T .A 



T" 

+ ( k s + nrf— 



(65) 



where 



A(r) = 2k s + nrf - 2nf . 



(66) 



Above equations show that two variable <fi and B express the perturbative variables and 
so these are master variables. Therefore, we must construct evolution equations for these 
variables in order to examine the stability of background solution. That for B has already 
been derived as f l5T|) . but this include metric perturbations. Then, substituting fl6"5"l) into 
fl57|) reads 



B - f 2 B" + f I In 



/ 



-5 



+ ^ — (ln(r n T'))' x (AnrEfB - 2nrfT$ - 2(k s T + nrfT')<p) = . (67) 



The evolution equation for is derived from nr f x (I62l) + ( l63l) with ( 165]) (see 
result is 

/V + / 2 ( In ' ^ 2 

/ 



0; 



and the 



nr 2 T 



r 2 fT 
(AT)' T" 
AT 



+ 



±K a {Er n ) | 2/ 2 (£r , 



(K s rT + nr 2 fT')-n(r 2 fT 

fr 



B = 



(68) 



nr n+1 T r n T \ V r "^ 2 (^ T ) 2 , 

These coupled two equations ( 167|) and (|68|) determine the behavior of and B. These 
two functions determine all perturbative variables, so these two equations are the master 
equations. 

Here, we derive a Schrodinger equation with two components like (I42p from these two 
equations. Against the case for vector type perturbations, it is more complicated because 
there is 0' in evolution equation for B. Therefore, in order to transform these equations into 
a Schrodinger type equation, we must eliminate this . For this, we must consider linear 
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combinations of the master variables </> and B. For example, the following combinations are 
fit for our purpose; 



rVT 7 



C 



1 



B 



n{r n E) 



A r ' ' r( n - 2 V 2 y/n(K a -n) V A 
By using these variables and tortoise coordinate r*, we can obtain 



n 



c 



where 



n 



V,(r) 7 c (r) 
VJr) V m (r 



and 



/ ( A( s r r" 



' ^4T T' ^ 



2n(r n EYf- 
ATr r 



fT' 



VJr) 



V em (r) 



k s — n \T' 
n r n l 2 A 

f r (n-2)/2Q r (f dt 
f 



r n ~ 2 (AT) 2 J J rVT 
A(r n E)f 2n(r n E)f 2 



AT ( rVT 



AT 



rT 



+ 



T 



In 



fT' 



1 



r r (n-2)/2 

A(Er n ) 2 



2n{r n E) 2 f 2 



In 



r n ~ 2 (AT) 
r 2n ~ 2 {AT) 2 



(69) 



(70) 



(71) 



(72) 



ATr^ 1 J r n AT \ \ f 

Here, we Fourier transforme like \& —> tye lU)t and £ C,e lu)t . Note that, in Einstein limit, 
we can decompose the above equation into two Schrodinger equations by taking the linear 
combinations like ip — a>±( where 



1 



a± 



A(r n E) 



n{n — 1) 



k, — n 



n + 1) \ M±\ M' 



16(k 8 -n){r n Ef 



n(n — l)(n + 1)' 



(73) 



This is consistent with Ref 



•B. 



We can do this because T(r) 



„n-l 



in Einstein theory. 



Against this, because T(r) is more complicated in general Lovelock theory, we must consider 
the above coupling system. 
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D. Condition for Instability 



In this subsection, we show the criterion for the instability under scalar type perturba- 
tions. 

Here, we show that "if 2T' 2 — TT" takes negative values somewhere in r > r out , charged 
Lovelock black holes have the instability" when T' is always positive. In order to show this, 
we here also define the inner product as 



dr* [v^vD 2 + ££ 2 ] 



(74) 



where ^> = (\1/, () T . For this proof, it is convenient to use the following inequality; for any 
test function ^test £ x C£°, the lower bound of spectra for % with x satisfies 



L0 2 ■ {$ test , $ test ) < {$ test , H^test) ■ 



(75) 



This inequality suggests that there exist instabilities if we can find trial function ^test which 
satisfies (^ tes t, WStest) < 0. 

We assume 2T' 2 — TT" has negative regions and define J as a closed set in the region 
2T' 2 —TT" < 0. Then, we choose a trial function as $test — (^o? 0) T where \l/o is a sufficiently 
smooth function with compact support on I. Using this test function, (^test, H^test) is 
evaluated as 



($ test , H% est ) = ^dr*[-%d 2 ^ + V g (r)\^ \ 2 } 



dr* 



d r *y + f[d r ln 



AT 
rVf~' 



f ( V T" 

+2. s - 2(. s -n)— --)|* 0| 



2n{r n E) 2 f 2 
ATr n 



r n ~ 2 (ATy 



(76) 



Note that we neglect the boundary terms in second equality because \l/ is zero at the 
boundary of /. Furthermore, using the relation 



AT = 2(k s - n)T - nr n+2 M'(r 
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the second line of the last equation can be evaluated as 



dr* 
dr* 
I dr* 



2k. 



f 



nr 



2(k„ - n) 



rpf rpf/ 



9 f 



AT 2T> 



2{k s - n) 



7 l*o| 



rpl rpf/ 



2k, 



nr y 2 (K s -n)-nr^^T 2T 
f ( T' T" 



|*r 



nr \T 2T' 



) l^ol 2 


= k s J dr* 







f 



nrTT 



_ (2T 12 -TT") |*o I 



(77) 



where we use the positivity of T and X" and also use eq. (IT8]) in the inequality. Then 
{^test, H^test) satisfies the following inequality; 



(Vtest, H^test) < / dr 



d r ^ + f[d r \n 



AT 

rVT' 



2n{r n Eff 2 



+k s / dr* 



ATr n 
f 



In 



*o 

fr 



r n ~ 2 (AT) 2 

— (2T /2 — TT") |*o 1 2 
nrTT' v 



l*r 



(78) 



In this equation, the second integral must be negative under our assumptions, therefore the 
second integral of (J78l) tends to — oo when k s = 1(1 + n — 1) — > oo. Against this results, the 
first integral of ( J78|) converge when I — > oo like 

2 



dr* 



d r *^o + f «9 r ln 



/ dr* 



AT 

rVT' 



d r ^ + f[d r hi 



T 



2n(r n E) 2 f 
ATr n 

2" 



hi 



fr 



r n ~ 2 (AT) 2 



l^ol 



rVT' 



(79) 



This is because integrands converge uniformly on I like 

AT 



<9 r »* + f \ d r hx 



rVT 1 



V ^d r .* + fld r hi 



T 



rVT' 



2n(r n E) 2 f 2 



In 



fr 



l*n| 2 ->■ 



(80) 



ATr n \~ \r n ~ 2 (ATf 

where these uniform convergences are supported by the theorem that continuity functions on 
closed set have maximum value and minimum value. Therefore, summarizing these results, 
we can denote that the r.h.s of ( 17H|) tends to — oo as £ — > oo. This means that the lower 
bound of the spectra is negative in sufficiently large i modes and so background solution 
has instability for these modes. 



25 



We'd like to summarize this section. We assume T' > in r > r out . Under this assump- 
tion, we derived master equations. We can unify these equations as a Schrodinger equation 
with two components. We show that this Schrodinger operator has negative spectra when 
2T' 2 — TT" has negative region. Therefore, we can denote that 2T /2 — TT" is crucial for the 
stability of charged Lovelock black holes. This criteria is same as neutral case 

We have not shown the inverse statement so far. Then, even if 2T' 2 — TT" is always 
positive, we can not declare this black hole is stable. For example, in Einstein case, T(r) is 
r n_1 , so 2T' 2 — TT" = n(n — l)r 2 ( n ~ 2 ) > 0. Then we can not say anything for Einstein case. 
Same as this, in 6-dimensions, this function can be evaluated as 



2 _ _ o r\p 2 M - 2r 5 ) 2 + 8foQ 2 (p 2 M + 3r 5 ) 

^ r9(l + 2/3 2 M(r)) >U 

so we cannot say anything for scalar perturbations in 6-dimensions . 



(81) 



VI. NUMERICAL RESULTS 

In section II III and section |V] we have shown that the behavior of T(r) is crucial for 
instability. In detail, T' is critical for ghost, T" is for tensor perturbations and 2T' 2 — TT" 
is for scalar perturbations. In this section, we check the behavior of these functions. For 
neutral cases, we can examine analytically because we can reduce these functions into the 



polynomial functions of ip [19j. However, such reduction can not be performed in charged 
case. Therefore, we numerically check the behaviors of T", T" and 2T' 2 — TT' for various 
(|Q|,.M) with some Lovelock couplings. 

In numerical calculation, we must use dimensionless parameters. So far we have discussed 
in the unit /3i = 1; in this unit, we can not fix the scale of length. Same as our previous 



paper 20|, we here also use (3 2 for fixing this scale. This constant has a dimension of 
length squared. Therefore, we can relate r, if), JH, Q and Lovelock couplings f3 m s to the 
dimensionless parameters as follows; 



r 2 

Hereafter, we show some results for 5 ~ 8 dimensions. The strategy of our numerical 



calculation is basically same as our previous paper except for checking 2T' 2 — TT" 20[ . Note 



that fi ex is the dimensionless extremal mass parameter which can be calculated from ( jTTl) . 
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^tensor is a border between stable and unstable for tensor type perturbations and ^scalar 
corresponds to that for scalar type perturbations. 

A. 5-dimensions 

As we have mentioned above, we use 02 for fixing the scale of length. Then we need not 
regard the Lovelock couplings in 5-dimensions. 

We present the numerical results for 5-dimensions in FigJH For this figure, we check the 
region where fi e x(\Q\) ~ ^ex(\Q\) + 3 for each \Q\ and the mesh size is dfi = dQ = 10~ 3 . 



As shown in our previous paper 19J, the tensor-unstable region lies thinly on the extremal 
mass fi ex (\Q\) and [itensor(\Q\) converges to 0.5 when \Q\ — )■ 0. In this limit, \i ex also converges 
to this value. Note that this thin region is over at \Q\ ~ 3 which has been checked in our 
previous analysis. Against the tensor-unstable region, for scalar type perturbations, the 
unstable mass range is relatively wide and this unstable region localizes near the /i-axis; the 
upper line of this region is approximately expressed as /i ~ 2.914 and this region suddenly 
disappear near \Q\ ~ 0.62. For the ghost instability, we have already checked the positivity 
of V in (|2"gj). 

Then, roughly speaking, charged Lovelock black hole with /i < 2.914 has the instability 
for scalar modes when \Q\ < 0.62; if 0.62 < \Q\ < 3, nearly extremal black hole is unstable 
for tensor type perturbations. We cannot detect the instability when black hole has more 
charge. 

B. 6-dimensions 

In 6-dimensions, same as 5-dimensional case, we need not alter Lovelock coefficients. 

As we have checked in section 1111} T' is always positive so there is no ghost. It has also 
mentioned in section [V] that IT' 2 — TT" is positive definite in 6-dimensions, which means 
we can not find the instability for scalar type perturbations. Then we can only detect the 



instability under tensor perturbations and the results are same as our previous analysis [20] ; 
same as the tensor-unstable region in 5-dimensions, this region exists just on fi ex (\Q\) and 
disappears at \Q\ ~ 3.28. However, there also exists a difference; ^tensor converges to 0.27 
in \Q\ — > while [i ex — > 0. Therefore, we can say that black hole with [i ~ /i ex has the 
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FIG. 4: Numerical results in 5-dimensions. The scalar-unstable region localizes near the //-axis and 
there exists a slight gap between ^tensor and fjL ex - This gap is about O(10 _1 ). Notice that "stable" 
means both T" and 2T' 2 — TT" have no negative region. This is same for the following figures. 

instability under tensor type perturbations when < \Q\ < 3.28. Especially, black hole has 
the instability also in neutral case against the 5-dimensional case. 




o £ 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 
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FIG. 5: Numerical results in 6-dimensions. Unlike 5-dimensional case, there are no scalar-unstable 
regions due to eq. (|8ip . The tensor- unstable region slightly lies on the extreme line /j, ex (\Q\). This 
region ends off at \Q\ ~ 3.28. Becasue fi ex — > as \Q\ — > 0, there also exists a tensor-unstable mass 
range in neutral case. 



28 



C. 7-dimensions 



We present the numerical results for 7-dimensions in Figj6j This figure is calculated with 
C3 = 0.2 in 7-dimensions. For this figure, we check the region where fJ* e x(\Q\) ~ HexilQl) +4.2 
for each \Q\ and the mesh size is d\x = dQ = 10~ 3 . 

This diagram is almost same as 5-dimensional cases; the tensor-unstable region clings to 
Hex{\Q\) and the scalar one does to //-axis. Our previous analysis shows that the tensor- 
unstable region exists in < \Q\ < 4.695 [20|. The upper bound of the scalar- unstable 
region is \i ~ 3.99 and this region vanishes at \Q\ ~ 0.516. When c 3 changes, the upper 
bound etc. change but the appearances of these unstable regions do not change. 

When C3 = 0.2, there are no ghost regions. However, when C3 is larger than 0.25, we can 
find ghost region near the origin of the diagram (see Fig. 14 of our previous paper 20] ). 

Therefore, same as 5-dimensions, we can roughly say that black hole suffers from the 
instability under scalar perturbations if Q ~ and \x is smaller than a certain value and has 
instability under tensor modes when \x is as small as fi ex . Against 5-dimensional case, there 
exists c 3 dependance for ghost regions. 



7D c 3 =0.2 near |Q|=0 
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FIG. 6: Numerical results for 7-dimensions. We calculate this with C3 = 0.2. Both ^.tensor and n ex 
converge to 0.667 as Q — > 0. Then there is no instability under tensor perturbations in neutral 
case. Note that there exists no ghost region in C2 = 0.2. The appearances of the diagrams do not 
strongly depend on C3 except for ghost regions. 
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D. 8-dimensions 



We present the two figures for 8-dimensions. FigJ7]is the numerical result for C3 = 0.07; we 
calculate the region in fi ex (\Q\) < ji < ^ e x{\Q\) + 0.05 and the mesh size is dfi = dQ = 10~ 4 . 
This figure is almost same as that for 6-dimensional case. In this figure, there are no 
parameters which make 2T' 2 — TT" negative, so we can only find the instability under 
tensor type perturbations. For tensor type perturbations, the unstable region slightly lies 
on n ex {\Q\) and there also exists a gap between ^tensor and [i ex in \Q\ — > 0. These properties 
are similar to 6-dimensional case. 

Against the above results, there exist a scalar-unstable region in FigJHJ This figure is 
calculated with c 3 = 1 and check the region in fi ex (\Q\) < ji < fi ex (\Q\) + 30 with the mesh 
size d/j, = dQ = 10 -3 . This figure is very similar to 5-dimensional diagram; a tensor-unstable 
region exists just on extreme line ix ex {\Q\) an d a scalar- unstable region localizes near /i-axis. 
Furthermore, there is no ghost region in C3 = 1. 

In FigJ7| and FiglHl there are no ghost regions. However, when C3 becomes larger than 
5.92, a ghost region appears near the origin of the diagram (see Fig.18 of 20j). Then, the 
appearance is similar to 7-dimensional case with sufficiently large C3. 

In 8-dimensions, as we have presented, the appearance of diagrams are very responsive 
to the Lovelock coupling c 3 . When c 3 is very small, this is same as 6-dimensional diagram. 
As C3 becomes larger, the looks of the diagrams change from 5-dimensional results to 7- 
dimensional diagrams with large C3. It is still an open issue why such dramatical changes 
occur in 8-dimensions. 



E. Summary of Numerical Results 

In this section, we have numerically checked the condition for the instability and ghost. 
Here, we summarize the results. 

We numerically examine the behavior of T(r) for various Lovelock coupling in 5, 6, 7 
and 8 dimensions and plot the results in Q — /i diagrams. From these results, we can read 
some common properties. The first is locations of unstable regions. The unstable regions 
for tensor type perturbations slightly lie on extreme line \x ex . The regions for scalar type 
perturbation, if exists, localize at the /x-axis and extend to the relatively large /1. The second 
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FIG. 7: Numerical results for 8-dimensions with C3 = 0.07. The appearance of this figure is very 
similar to the result for 6-dimensions; There only exists the tensor-unstable region near extreme 
line /ji ex {\Q\). ^tensor converges to 0.0055 and fi ex tends to when \Q\ — > 0. Then there also exists 
the instability under tensor type perturbations when black holes are neutral. We cannot find ghost 
regions and scalar-unstable regions in C3 = 0.07. 




FIG. 8: Numerical results for 8-dimensions with C3 = 1. Against C3 = 0.07 case, there exist a 
scalar-unstable region near \Q\ = 0. Then, in 8-dimensions, which type instability black holes 
suffer from is very sensitive to third order Lovelock coupling C3. In this figure, \x ex — > and 
^tensor — >• 0.459 as \Q\ — > 0. Note that we cannot find ghost regions in C3 = 1. 
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is manner of tensor-unstable region near \Q\ = when there is no scalar-unstable region. 
For this case, in our numerical calculation, there must exist slight mass range in which black 
holes have the instability under tensor perturbations. Therefore, whether scalar-unstable 
region exists or not, when black holes with nearly extreme mass have slight charge, they 
must be unstable; which type instability they have depends on parameters and dimensions, 
but they have at least one type instability. 

These results also lead some open questions. The first is dimensionality. The diagram 
for 5-dimensional case and that for 7-dimensional case are very similar: the result for 6- 
dimensions and that for 8-dimension with small C3 are alike. These remind us that the 
behaviors of black holes against perturbations are different in even dimensions and odd 
dimensions. However, we do not have the answer for such a dimensionality. The second is the 
response to variation of Lovelock couplings. In 8-dimensions, for example, the appearance 
of diagram is sensitive to third order coupling C3. These are still open questions whether 
this is proper to 8 dimension or not and what causes such peculiarity. 

VII. CONCLUSION 

We have studied the stability of charged Lovelock black hole solutions. We have derived 
master equation for vector type perturbations and scalar type perturbations. These are 
the Schrodinger equations with two components. For vector type perturbations, we have 
shown that the Schrodinger operator for this type perturbations is an essentially positive 
definite self adjoint operator. Then charged Lovelock black holes are stable under vector type 
perturbations. On scalar type perturbations, we have presented the condition for instability. 
In detail, if IT 12 — TT" has negative regions, charged Lovelock black holes are unstable under 
scalar type perturbations. For tensor type perturbations, we have already shown that T' 
and T" are crucial for stability. By numerically checking these criteria, for example in the 
2nd order and the 3rd order Lovelock theory, we have shown that there exists an unstable 
parameters; nearly extremal black holes have the instability under tensor type perturbations, 
and black holes with small charge have the instability under scalar type perturbations even 
if black holes have relatively large mass. 

One of future works is the exploration of more general conditions for stability under scalar 
type perturbations. In this paper, we have shown that black holes are unstable if 2T' 2 — TT" 
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has negative region. However, the inverse statement has not been proved. Hence, so far, we 
can not say anything when 2T' 2 — TT" is always positive. Furthermore, by this criterion, 
we can not detect the instability under scalar type perturbation in Einstein theory 9[. In 
this sense, it is interesting to find out more general conditions. 

It is interesting to investigate the relation between dynamical instability we have shown 
in this paper and thermodynamics. On the thermodynamics for Lovelock black holes, varia- 



231 ] . In our paper, we have found that 



tion of Lovelock coefficients is also examined in Ref. 
the appearances of Q — \i diagrams change dramatically as Lovelock coupling C3 varies in 
8 dimensions. Then it is interesting if such dramatical changes are found in the thermody- 
namics. 

The relation between instability and gravitational collapses might be important. In Love- 
lock theory, collapses of dust clouds have been examined 24|. Furthermore, these are ex- 



tended to charged dust clouds [25|. In these, dependence of dimensions are found for, for 
example, naked singularity formations. Our results also depends dimensions, so there may 
exit relations between the instability and the gravitational collapse in Lovelock gravity. In 
dust collapses, the authors of above papers have also pointed out the tendency that higher 
curvature collections suppress formations of apparent horizons. These results should express 
that higher curvature collections make attractive force weaker, and this property might be 
related to the instability of Lovelock black holes we discussed in this paper. 
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